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We study the Josephson junction between two rf-wave superconductors, which is discussed as an imple- 
mentation of a qubit. We propose an approach that allows to calculate the decoherence time due to an 
intrinsic dissipative process; quantum tunneling between the two minima of the double-well potential excites 
nodal quasiparticles which lead to incoherent damping of quantum oscillations. The decoherence is weakest in 
the mirror junction, where the contribution of nodal quasiparticles corresponds to the superohmic dissipation 
and becomes small at small tunnel splitting of the energy level in the double-well potential. For available 
experimental data, we estimate the quality factor. 

PACS: 85.25.Cp, 85.25.Hv, 74.50.+r, 73.23.-b 



Among various candidates for physical implemen- 
tation of quantum bits, solid-state proposals, and in 
particular superconducting devices, have a number of 
advantages, e.g., scalability and variability p. Particu- 
larly interesting are the so-called quiet qubits, which are 
intrinsically degenerate, i.e., do not require any exter- 
nal source for maintaining the degeneracy. Such qubits 
can be realized in systems involving d-wave supercon- 
ductors [2]- Recently, it was experimentally demon- 
strated that a double-well potential is indeed realized 
in the Josephson junctions between d-wave supercon- 
ductors The qubit variable in this case is the phase 
difference ip across the junction. The energy of the phase 
qubit has two nontrivial minima as a function of the 
phase difference (see Fig^. Alternatively, a quiet flux 
qubit can be realized if the spontaneous flux is gener- 
ated in the loop of d-wave superconductors 0] . The two 
qubit implementations are quite similar; for definiteness 
we shall speak about the phase qubit. 

In such intrinsic qubits, there are also intrinsic mech- 
anisms of decoherence even at low temperatures. The 
quantum tunneling of the phase between the two min- 
ima leads to fluctuating voltage across the junction, 
which excites quasiparticles. The dissipative current 
across the interface arises, leading to a finite decoher- 
ence time Tip. The knowledge of is essential for es- 
timating the efficiency of the qubit: short decoherence 
time makes the qubit senseless, while a long enough de- 
coherence time opens the way for quantum correction al- 
gorithms that in principle allow to perform an infinitely 
long computation 
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FigQ Schematic dependence of the Josephson energy 
U on the phase difference (p (in the flux qubit ip is sub- 
stituted by 27r<l?/<I>o, with Oo the flux quantum). The 
barrier of the height 2Ej separates two nontrivial min- 
ima. The splitting of the lowest energy level due to the 
tunneling across the barrier is denoted At. 

The relevance of the quasiparticle processes at low 
temperatures is specific for d-wave superconductors. In 
the conventional s-wave case, the quasiparticle trans- 
port below the gap is suppressed. At the same time, 
in gapless anisotropic superconductors the gap vanishes 
in certain directions (the nodal directions), hence the 
low-energy quasiparticle appear. In the present letter, 
we consider a DID Josephson junction (D = rf-wave su- 
perconductor, I — insulator), and study the decoherence 
due to nodal quasiparticles (quasiparticles moving along 
the nodal directions). 

Decoherence time (general strategy). Theo- 
retical description of the quantum dynamics of a tun- 
nel junction between two s-wave superconductors was 
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developed in Ref. IBj (see Ref. j7j for a review). The 
effective action for the phase difference if was obtained. 
Later this description was generahzed to the case of d- 
wave superconductors in Refs. |H1I1|. The effective ac- 
tion for if is similar to the general case considered by 
Caldeira and Leggett | lUllllj . who studied influence of 
dissipation on quantum tunneling in macroscopic sys- 
tems. The dissipation was described as being due to the 
interaction with a bath of oscillators (the environment) . 
The "strength" of the environment, depending on the 
frequency lo, is characterized by the spectral function 
J{uj). In the Josephson junction, the environment is 
represented by the quasiparticles, and the spectral func- 
tion is given by hlihuj / e) / e, where / is the dissipative 
quasiparticle current taken at "voltage" Huj/e 6 . 

A system living in a double-well potential and de- 
scribed by an extended coordinate can be "truncated" 
to the two-state system (spin 1/2) with the two states 
= ±1) corresponding to the minima of the poten- 
tial (see Fig^l. The theory of dissipative two-state sys- 
tems is thoroughly elaborated jjl] for the cases when 
the spectral function behaves as J(w) oc oj^ up to some 
high-frequency cutoff. The situations when s — 1, s > 1, 
and < s < 1 are called ohmic, superohmic, and sub- 
ohmic, respectively. In this language, the dissipation 
due to nodal quasiparticles in the Josephson junction is 
superohmic, as we demonstrate below. 

What is the decoherence in such a system? Assume 
that during the time t < the system is held in the right 
well (i.e., at az = !)• At i = the constraint is released, 
and we consider the expectation value of the system co- 
ordinate: P{t) — {az{t)). Below we shall encounter the 
superohmic case at zero temperature. Then 



Pit) = cos(Att/;i) exp(-i/r<p) 



(1) 



— the cosine describes coherent oscillations between the 
two wells (At is the tunnel splitting of levels, see Fig^J 
while the exponential leads to their incoherent damping. 

The decoherence time T^p is expressed in terms of the 
spectral function Returning from the general the- 
ory to the particular case of the Josephson junction, we 
write the corresponding result as 

_ 4e _ AttH 

" V^(At/e) " Vei?,/(At/e)' ^ ' 

where Sip is the distance between the potential minima 
and Rq — hjle^ w 13kf2 is the quantum resistance. 
Comparing the decoherence time with the characteristic 
time of oscillations between the wells, H/ l^t, we obtain 
the quality factor 



Q 



27rAt 



2h ^ S(p^eRgI{At/e)' 



which must be large for successful operating of the 
qubit. 

In the DID junction, the tunnel splitting A^ is much 
smaller than the order parameter A, hence is deter- 
mined by the quasiparticle current at low "voltage" . 

Quasiparticle current. Motivated by the ex- 
periment P], we consider the grain-boundary Joseph- 
son junction between two quasi-two-dimensional d^^-y^- 
wave superconductors with cylindrical Fermi surfaces. 
The orientations of the superconductors are character- 
ized by the angles between the a-axes and the normal to 
the interface (the x-axis) — see Fig|2 According to Ref. 
|3|, we consider the mirror junction, in which the mis- 
orientation angles on both sides are equal in magnitude 
but opposite in sign, a/ — a (we take —45° ^ a ^ 45° 
because all physically different situations in the mirror 
junction are realized in this interval). The order pa- 
rameter depends on the direction (parametrized by the 
angle 9) and the distance to the interface: 

AlM^, 9) = AlA^V^^'" cos {2{e T a)) , (4) 

where the indices L and R refer to the left- and right- 
hand side of the junction, respectively. 

The quasiparticle current in the tunneling limit at 
low temperatures, fcsT <C is given by 



I{Huj/e) 



1 r\,m^r,^^ 



sRn J-tt/2 D jo 

xN{E-huj,9)N{E,9) 
d9D{9) cos 9. 

J -it/2 



(5) 



Here R^ is the normal-state resistance of the interface, 
N{E, 9) is the density of states (DoS) at the interface, 
normalized to the normal-metal value, and D{9) is the 
angle-dependent transparency of the interface. We have 
not labelled the DoS by the indices L and R because 
Nl{E,9) — Nii{E,9) in the mirror junction. 

Below we calculate the nodal contribution to the cur- 
rent Q at fiij <C Aq, where Aq — A(±oo) is the bulk 
amplitude of the order parameter. The angle integration 
contributing to the current is then limited to narrow an- 
gles around the nodal directions, where the low-energy 
DoS is nonzero (as we shall see below, the width of the 
angles is S9 = Huj/Aq). 

To calculate the DoS, we employ the quasiclassical 
approach. The quasiclassical matrix Green function 
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FiglH DID junction of mirror orientation a/ ~ a. An 
electron e moving along a truly nodal direction 6i of the 
left superconductor, tunnels into an induced nodal di- 
rection of the right superconductor. A[i{9i) ^ 0, there- 
fore the electron experiences the Andreev reflection; the 
hole h returns to the interface, and after reflection at 
the interface escapes into the bulk along the truly nodal 
direction —Oi. In this process, the total current into 
the bulk of the right superconductor is composed of the 
Cooper pair along 9i and the hole along ~6i. 

obeys the Eilenberger equation |21 and satisfies the nor- 
malization condition = 1. It can be parametrized as 
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1-ab 
l + a6' 
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2a 



l + ab' 
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2b 



1 + ab' 



(7) 



then the normalization condition is automatically sat- 
isfied. The equations for the new functions a{x, 9) and 
b{x, 6) take the form of the Riccati equations llBj : 

hvF cos9{da/dx) - 2iEa + A*a^ - A = 0, 

hvF cos e{db/dx) + 2iEb - Ab^ + A* 0, (8) 

where vp is the absolute value of the Fermi velocity vp, 
and 9 denotes the angle between wp and the x-axis. 

In the tunneling limit, the DoS is calculated at 
an impenetrable interface. Let us consider, e.g., the 
right superconductor (the right half-space). We need 
to find the low-energy DoS in two cases: 1) in the 
vicinity of a nodal direction, so that E,A{9) <C Aq, 
2) at a gapped direction, so that E ^ A{9). In 
the first case, the space scale £,e = hvp cos9/\£+\ on 
which the quasiclassical Green functions vary (we de- 
note £± — \J E"^ — |A(oo, ±6')!^), is much larger than 
the coherence length ^ = hvp/2TrkBTc on which varia- 
tions of A occur. This allows us to regard A as constant 
when integrating Eqs. (O over x. In other words, the 
functions a and b at low energies do not feel the suppres- 
sion of A near the interface, because it takes place on a 
small scale. In the second case, the spatially dependent 



parts of a and b are proportional to E/A{9) ^ 1 and 
hence small. Thus a and 6 at the interface are equal to 
their bulk values, as if A was constant. 

Thus we can regard A{x,9) as equal to the bulk 
value Ao{9) = A{oo,9). The integration of the func- 
tions a and b over x in Eqs. ((HJ is stable only in the 
directions determined by the sign of cos 0. At cos 9 > 0, 
the function b(x, 9) is stably integrated from x = cxd to 
the interface {x — 0), hence 

b{0,9) ^ b{oo,9) = i{E ~£+sgnE)/Aoi9). (9) 

At the same time at cos 9 > 0, the function a is stably 
integrated from the interface to a; = oo. Therefore to 
find a{0,9), we consider the trajectory directed along 
TT — 9. Since cos(7r ~ 9) < 0, the function a is stably 
integrated from a; = oo to the interface. Finally, the 
direction tt — is converted to 9 upon reflection at the 
specular interface: 



a(0, 9) = a(0, tt - 9) = a(oo, tt - 9) = 
= i{E-£^ sgnE)/A*{-9). 

As a result, the DoS N — Keg at the interface is 

\E\ {£++£_) 



N{E,9) = Re 



E^ - Ao{9)Al{-9)+£+£^ 



(10) 



(11) 
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The gap in the spectrum 
min(|Ao(0)|,|Ao(-0)|). 

The DoS is symmetric, N{9) — N{~-9), because the 
Green functions are continuous upon reflection. Thus in 
each superconductor there are two "truly" nodal direc- 
tions 9i (i = 1,2) in the interval —tt/2 < 9 < 7r/2, and 
also two "induced" nodal directions — 6*.;. Near a nodal 
direction Eg{9) = 2Aq\9 — 9i\. Along a truly nodal di- 
rection, the gap vanishes and the DoS is the same as in 
the normal metal, N{E) — 1. For an "induced" nodal 
direction this is so only near the interface. 

In the left superconductor, the truly nodal directions 
are 9i^2 = a ± 45°. Due to the mirror symmetry, the 
truly nodal directions of the right superconductor coin- 
cide with the induced nodal directions of the left one, 
and vice versa. In total, there are four nodal directions 
in the junction, which are symmetric with respect to the 
interface normal. 

In this situation, the transport is due to the pro- 
cesses of the following type. An electron moving along 
a truly nodal direction 9i of the left superconductor, 
tunnels into an induced nodal direction of the right su- 
perconductor (see FiglSJ. However, the electron cannot 
escape into the bulk of the right superconductor because 
Aii{9i) ^ 0. Therefore the electron experiences the An- 
dreev reflection; the hole returns to the interface, and 
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after reflection at the interface escapes into the bulk 
along the truly nodal direction —9i. In this process, the 
total current into the bulk of the right superconductor 
is composed of the Cooper pair along 9i and the hole 
along —9i, which is overall equivalent to the transfer of 
one electron. 

The nodal contribution to the current (O appears 
only due to integrating in the vicinity of the nodal di- 
rections where Eg < hio. The DoS near the nodal di- 
rections at small energies can be found from Eq. 
Below we distinguish the general case when Ao(0) ^ 
±Ao(— 6*), and two special cases: Ao(6') = Ao(— 0) (at 
a = 0°) and Aoie*) = -Ao(-6l) (at a = 45°). 

At a = 0°, the truly nodal and induced nodal di- 
rections coincide in each superconductor, and Eq. 
yields the BCS-like DoS: 



Noo{E,e) = Re(|i?|/^£;2 - |Ao(0)|'). (12) 

At a = 45°, the truly nodal and induced nodal di- 
rections again coincide, and Eq. Hll|l yields the DoS of 
the inverse BCS type: 



iV45o(i?,e) -Re(7i?2 - \^o{ef/\E\). (13) 

Finally, if |a| > /iw/Aq and 45° - |a| > ftw/An (i.e., 
a is not too close to 0° and ±45°), then Ao(6') in the 
essential angle of the width 59 — huj/Ao around a nodal 
direction is much smaller than Aq{—9). Then in the 
region of energies that contribute to the quasiparticle 
current, |Ao(6')| < \E\ < Hlu <€. |Ao(-6')|, the DoS is 
again given by the inverse BCS formula: 



Ng{E,9^ 9,) = Re - |Ao(^)l7l-B 



(14) 



Figure |31 demonstrates the DoS at different an- 
gles 9, which are parametrized by different ratios 
Aq{~9)/Aq{9). At Aq{-9) ^ Ao{9), the DoS has the 
BCS-like square- root singularity near Eg [see Eq. (fT^ ]. 
At Ao{-9) ^ Ao{9), the DoS has the inverse-BCS be- 
havior near Eg [see Eq. (I14f) ]. 

Inserting Eqs. l(T^ - ((n|l into Eq. ©, we obtain: 



I{huj/e) 



cRn Ao 



E 



D{9,)cos9, 

b 



(15) 



where 9i 2 = a ± 45° and A is a number, which de- 
pends on the orientation of crystals: A{Q°) w 0.46, 
A(45°) « 0.19, and A{a) = 2A(45°) « 0.37 when a 
is not too close to 0° or ±45°. 

In Refs. miSl, the quadratic current -voltage charac- 



teristic, / 



was obtained for the case of aligned 
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FigEl Density of states following from Eq. (IIH . The 
energy is normalized to Ao(S), while Ao(— ^) is varied. 



Estimate. Equations Q, H15() yield: 



27r 



^5ip^A{a) D{9,)cos9i/D^ 

1=1.2 



To proceed further, we need to estimate the tunnel 
splitting At (see Fignj. For the estimate, we assume 
that the second harmonic dominates in the energy- 
phase relation, U{(p) = Ej{l + cos2(^), and the energy 
of the levels is small compared to Ej . Then the tunnel- 
ing action is calculated between the points tp = —tt/2 
and 7r/2, and we obtain: 



At = (4^2EjEc/tt) exp(~^/2Ej/E, 



nodal directions (i.e., for the a/ a orientation). 



(17) 



where Ec = e'^/2C is the charging energy (C is the 
capacitance of the junction). 

To obtain a numerical estimate, we take the charac- 
teristics of the junction as in the experiment of Il'ichev 
et al. inj. The capacitance of the junction is C ~ 10~^^F 
|14j . hence Ec/ks ~ O.IK. The characteristic Joseph- 
son energy is of the order of several Kelvin. For es- 
timate, we take 2Ej/kB = 7K. The resistance of the 
interface is Rn ~ 50 O (2] . 

As a result, At/Zcs ^ 2.5 • 10^"* K. Finally, we esti- 
mate 5lp ~ TT, Ao/ks ^ 200 K, and assume a thin 6- 
functional barrier with transparency D{9) = Dq cos^ 9, 
then the quality factor is Q ~ lO'^ ^ 10"'. Here we have 
retained only the order of magnitude for Q, because we 
cannot expect a higher accuracy in the case when im- 
portant characteristics of the junction (e.g., C and Ej) 
are known only by the order of magnitude. We also 
made an essential assumption that the second Joseph- 
son harmonic dominates. 
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The latter assumption can be realized under spe- 
cial conditions, while in a more common situation the 
first and the second harmonics are of the same order. 
Estimates for this case were done in a recent work 
[T5) . where the characteristics of mesoscopic junctions 
between high-Tc superconductors were experimentally 
studied and theoretically analyzed. A characteristic 
value of At ~ O.IK was reported under the conditions 
that correspond to Rn ^ 100 ft. Assuming such param- 
eters for the mirror junction, we obtain Q ^ 10 10^. 

The above estimates for Q are very different. At the 
same time, a general consequence of Eq. H16|l is that the 
quality factor grows as the splitting At becomes smaller. 
We note in this respect, that the values of the critical 
current (and hence the Josephson energy) measured in 
Refs. IHIE], are much smaller than expected. If the 
critical current is enhanced to the expected value, then 
At decreases, which finally leads to an increase of Q. 

If a / 0°, the low-energy quasiparticles are pre- 
sented not only by the nodal quasiparticles, but also 
by the midgap states (MGS) with zero energy IBj. In 
the case of specular interface and clean superconduc- 
tors, considered in this paper, the DoS corresponding 
to the MGS is proportional to S{E), hence the MGS on 
the two sides of the interface do not overlap and do not 
contribute to the current at a finite voltage. 

In the asymmetric case, when aL ^ ±afl (precisely 
speaking, when ||aL| — |a/j|| > huj/Ao), the nodal direc- 
tions of the left and right superconductors do not match 
each other. Then the transport from nodal to nodal di- 
rection is suppressed. However, a more important trans- 
port "channel" arises: between the nodal directions and 
the MGS. This leads to a stronger decoherence than in 
the symmetric case. 

In the mirror junction, the MGS contribute to the 
quasiparticle current if they are split and/or broad- 
ened To take into account the contribution of the 
MGS into decoherence, the present approach should be 
considerably modified. This issue requires a separate 
study. 

In conclusion, we have proposed an approach that 
allows to calculate the decoherence time due to nodal 
quasiparticles in the DID junctions, which can be used 
as phase or flux qubits. The dissipation in the mir- 
ror junctions is weaker than in the asymmetric ones. 
We find the superohmic dissipation with s = 2 in the 
mirror junction, which becomes weak at small tunnel 
splitting of the energy level in the double-well poten- 
tial. For available experimental data, we estimate the 
quality factor. The superohmic case is most favorable 
(compared to ohmic and subohmic) for possible qubit 
applications. 
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